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Abstract 

The quantum-resonant evolution of the mean kinetic energy (MKE) of the 
kicked particle is studied in detail on different time scales for general kicking 
potentials. It is shown that the asymptotic time behavior of a wave-packet 
MKE is typically a linear growth with bounded fluctuations having a sim- 
ple number-theoretical origin. For a large class of wave packets, the MKE is 
shown to be exactly the superposition of its asymptotic behavior and transient 
logarithmic corrections. Both fluctuations and transients can be significant 
for not too large times but they may vanish identically under some conditions. 
In the case of incoherent mixtures of plane waves, it is shown that the MKE 
never exhibits asymptotic fluctuations but transients usually occur. 
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The recent extensive studies of quantum resonance (QR) and related phenomena exhib- 
ited by the periodically kicked particle (KP), either in the presence or in the absence of 
gravity, have led to interesting new connections between the classical and quantum dynam- 
ics of nonintegrable systems [1-4]. These phenomena are now recognized as highly efficient 
methods for imparting large changes in the kinetic energy of cold atoms in atom-optics ex- 
periments [1,2]. The quantum KP in the absence of gravity is described by the Hamiltonian 



where (x, p) are position and momentum operators, k is the nonintegrability parameter, 
V(x) is a periodic potential, t' and t are the continuous and "integer" times, and r is the 
kicking period; the units are chosen so that the particle mass is 1, h — 1, and the period of 
V(x) is 2n. The most investigated case of QR, which is also the easiest to observe experi- 
mentally, corresponds to integer values of r/(2n). It was shown [3-5] that in this case QR 
manifests itself in an asymptotic linear growth, as a function of time t, of the mean kinetic 
energy (MKE) of either a general KP wave packet or an incoherent mixture of plane waves. 

Experimental observations of the QR phenomenon are, of course, limited in time. In ad- 
dition, since integer values of r/(2n) cannot be precisely realized experimentally, the growth 
of the MKE may stop or saturate after some finite time. This was shown in Ref. [4] when 
studying the evolution of an incoherent mixture of plane waves for r = 2n + e (|e| <C 1), 
using the standard potential V(x) = cos(x). It was found that the MKE (E t , e ) of this mix- 
ture essentially coincides with that for e = in the time interval < t < \ke\~ l l 2 , where 
(E ti e) ~ (Etfi) is already approximately given by the asymptotic linear behavior of (E tt0 ); 
however, for t > |A;e| -1 / 2 , the growth of (E t:t ) is suppressed due to dynamical localization. 
In Ref. [5], the momentum probability distribution of a wave packet for a two-harmonic 
potential was shown to be also robust, on some initial time interval, under small deviations 
e of r around r = 2n. These facts motivate the investigation of quantum-resonant evolution 
on arbitrary time scales, not just on asymptotic ones. In particular, since the asymptotic 
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behavior may start in practice at times t not too long, it is important to inquire about 
the exact nature of this behavior, i.e., whether fluctuations take place around the (average) 
linear growth for general potentials V(x) and what is the strength of these fluctuations. 

In this paper, the quantum-resonant evolution of the KP is studied in detail on differ- 
ent time scales for integer t/(2tt) and general V(x). In the case of the MKE (E) t of wave 
packets, we show that asymptotic fluctuations of (E) t occur only if V(x) is multi-harmonic 
and the initial wave packet is nonuniform (in a well-defined sense). These fluctuations and 
the average linear growth have essentially the same origin, of a simple number-theoretical 
nature. The fluctuations have bounded variation in time and are therefore much smaller 
than (E) t for large t. For a large class of initial wave packets, we derive exact closed ex- 
pressions for (E) t : (E) t is just the superposition of its asymptotic behavior (linear growth 
and fluctuations) and transient logarithmic [ln(i)] corrections. The fluctuations are exactly 
periodic or quasiperiodic in time and both they and the transients can be significant for 
not too large times. The transients vanish if V(x) and the initial wave packet satisfy some 
symmetry conditions. Then, the asymptotic behavior starts at t — 1. If also the fluctuations 
vanish, e.g., for mono-harmonic V(x), (E) t is an exactly linear function of t. In the case 
of incoherent mixtures of plane waves, we show that the MKE never exhibits asymptotic 
fluctuations but transients occur for nonuniform mixtures. 

Background. We first summarize known facts [3-5] about the KP quantum dynamics. 
The one-period evolution operator for (1), from t! = t + to t' = t + r + 0, is given by 

U = exp [-ikV(x)) exp (-irp 2 /2) (2) 

(H = 1). Because of the 27r-periodicity of (2) in x, it is natural to study the quantum 
evolution in terms of Bloch functions (pp(x) = ex${i(3x)il)p(x), where (3 is a quasimomentum, 
< (3 < 1, and if)p{x + 2ir) = ipp(x). In fact, the application of U on ip^x) results in a 
Bloch function associated with the same quasimomentum, U(pp(x) = exp(i{3x)i(}p(x); here 
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4>'p{x) is the 27r-periodic function ip'^x) = Up^p{x), where 

Up = exp [-ikV(x)\ exp [-it (p + (3) 2 /2\ . (3) 

Since the operator (3) acts only on 27r-periodic functions il)p(x), one can interpret x as an 
angle and p as an angular-momentum operator with integer eigenvalues n. Then, (3) may 
be viewed as the one-period evolution operator for a "/9-kicked rotor" (/3-KR). Given an 
arbitrary KP wave packet ty(x) and using the relation between ty(x) and its momentum 
representation ^f(p), it is easy to see that ty(x) can be always expressed as a superposition 
of Bloch functions, ^(x) = Jq dpexp(if3x)tf)p(x) : where 

ipffjx) = j — ty(n + (3) exp(inx). (4) 

\J2ll n 

One then gets the basic relation 

* t (x) = U^(x) = £ dPexp(i(3x)U{,ipp(x) (5) 

for integer time t, connecting the quantum dynamics of the KP with that of (3-KRs. The 
evolution of a /3-KR wave packet ipp(x) under the QR condition t/(2tt) = I, where I is a 
positive integer, was studied in Ref. [5] for general 27r-periodic potential V(x), 

V(x) = V m exp(-imx). (6) 

m 

It was found that, up to a nonrelevant constant phase factor, one has [5] 

ip/3,t( x ) = u piP/3( x ) = ex P [-ikVi3, t (x)\ ipp(x - trp), (7) 
where Tp = 7rZ(2/3 + 1) and 

AU ' h m sin(mr /3 /2) 1 ; 

General exact expressions. The expectation value (E) t of the kinetic energy p 2 /2 in the 
wave packet (5) can be expressed in the p representation with p decomposed into its integer 
(n) and fractional (f3) parts, p = n + (3. The resulting expression can be compactly written 
in terms of the Bloch functions ipp tt (x) = exp(i(3x)ipp it (x), using (4) (see note [7]): 



(E) t = U 1 dPj2( n + P) 2 \Mn + P)\ 2 = U L dp [ 
2 Jo „ 2 Jo Jo 



- d/3 dx 
2 Jo Jo 



dyp,t{ x ) 



dx 



(9) 



Under QR conditions, rj (2ir) = I, we see that Eq. (7) is satisfied with ipp replacing ipp on 
both sides, up to the factor exp(if3trp). This factor, however, will not appear in (9). We 
then easily find that (E) t is the sum of three terms, (E) t = (£) t (1) + (E)f ] + (E)f \ where 



k 2 r 1 r 2 ^ 



(pp(x)- 



dVp/x + tr, 



dx 



{E)f ] = k C dp r dxlm 
Jo Jo 



ipp(x) 



dx 



dx 



(10) 



(11) 



d<Pp( 



X 



dx 



(12) 



The constant term (12) is just the initial value (E) of (E) t . We now derive more explicit ex- 
pressions for (10) and (11). We then show that (10) is totally responsible for the asymptotic 
behavior of {E) t , including characteristic fluctuations. The term (11) is shown to contribute 
transient effects which can be significant in some cases. 



To express (10) more explicitly, we use an expansion following from Eq. (4): 

= IV^) | 2 = ^-^Cp{m) exp(imx), 
2ix m 



(13) 



where Cp{m) are correlations of the KP wave packet in angular-momentum (fixed (3) space, 



Cp( m ) = *( m + n + P)^*( n + P) = Y1 Cj(m) exp(2irijj3). 



(14) 



The last expression in (14) is the Fourier expansion of C 13(171), based on the obvious period- 
icity Cp + i(m) = Cp(m). We insert Eqs. (13), (14), and (8) in (10) and use the identity 



sm(mrpt/2) 



t-i 



cxp[im(t — l)rg/2] exp ( —imsTp ) . 



(15) 



sm(mrp/2) s=0 

Since Tp = nl{2j3 + 1), the integrations in (10) can be explicitly performed and we finally 
obtain the general exact expression 
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E mm ' W E (-l) ?(mS mS) G(ms-m' S ')( m - m/ )- 



(16) 



m,m' s,s'=l 

Consider now the term (11). Using (pp !t (x) — ex P(^^)'0/3,t( a:; ) with (4), we find that 



Im 



= -— E *(« + + /9)(n + ra' + 2/3) exp[i(n - n')4 (17) 



Let us introduce the quantity 



Gp{m) = E + n + PW(n + (3)(n + (3) =J2 G j( m ) exp(2mj (3) , 



(18) 



where the last expression is the Fourier expansion of Gp(m), based on the periodicity 
Gfs +1 (m) = Gp{m). Using Eqs. (17), (8), (15), (14), and (18) in (11), we get, after a 
straightforward calculation, 



(E)f = -i E V m E (-!)' mS Wc lms {m) + 2mG lms (m) 



(19) 



" m s =l 

The growth of (19) with t is always slower than linear since both Cj(m) and Gj(m) decay 
with \ j\. For sufficiently smooth ^f(p), this decay is fast enough to yield a finite value of (19) 
for all t. Also, (E)^ may vanish identically if some symmetry conditions are satisfied. Let 
have a definite parity, = ±§(p). Then, from (14) and (18), Cp(m) = C^(-m), 

i.e., Cj(m) = C-j(-m), and Gp{m) = —G-p(—m), i.e., Gj(m) = —G-j(-m). If also V(x) 
is odd, V(—x) = —V(x), i.e., U_ m = — V m , it follows immediately that {E)f^ = in (19). 



Asymptotic time behavior. We now show that the asymptotic behavior of (16) for large 
t is typically a linear growth with bounded fluctuations; then, since (19) grows slower than 
linearly (see above), the asymptotic average growth of (E) t is linear. Let us first change 
variables in the second sum in (16) from (s, s') to integer variables (a, b) defined as follows. 
We denote by g = g(m, m') the greatest common factor of (\m\, \m'\) and write m = gm , 
m! = gm' Q , where (m , m' ) are coprime integers; thus, ms — m's' = g(m s — m' s') in 
(16). For any integer a, a solution (s, s') of the Diophantine equation m s — m' s' = a is 
s = as and s' = as' , where (s , s' Q ) are coprime integers satisfying m s — rn' s' = 1; such 
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integers always exist for coprime (m , m' ), see, e.g., Ref. [6]. The general solution (s, s') 
of m s — m' s' = a is then 

s = as Q + bm' , s' = as' Q + bm (20) 

for all integers b. By the change of variables (20), l(ms — m's') = gla in (16) and, due to 
the decay of Cj(m) with \j\, C g i a (m — m!) is negligible if \a\ is too large. Now, for given a, 
the values of b in (20) are restricted by the condition 1 < s, s' < t in (16). If (m , m' ) have 
the same sign, the number JVj, of 6-values for t ^> max(|aso|, |asol) ^ s approximately [t/m], 
where m = max(|mo|, \m' Q \) and [x] denotes the integer part of x. The case of m and m' 
having different signs can be ignored since JVj, is much smaller than [t/m] in this case. One 
can therefore approximate (16) for large enough t as follows: 



Z mm'>0 a 1771 



(21) 



We decompose [t/m] as [t/m] = t/m + F t (m), where F t (m) = [t/m]— t/m are "fluctuations" 
of a simple number-theoretical nature. Correspondingly, (E)^ in (21) will be decomposed 
into a linear-growth part {E)f' L) and a fluctuating part {E)f' F) : (E)f ] w (E)f' L ^ + (E)^' F \ 
The expression for {E)f' L \ given by Eq. (21) with [t/m] replaced by t/m, can be easily 
shown to be precisely the known one for the asymptotic linear growth derived in Ref. [5] 
[see Eqs. (17) and (19) there]. In fact, from Eq. (14) we find by simple algebra that 

J2(-l) 9la C gla (m-m>) = ^Y;C^ g (m-m'), (22) 

a 9 L r =0 

where (3 r>g = r/ (gl) — 1/2 mod(l), r = 0, . . . , gl — 1, are resonant values of (3 [5]. Using (22) 
and the identity 2mm' /(gm) = \m\ + \m'\ — \m — m'\, (E)^ 1 '^ reduces to the expression (17) 
with (19) in Ref. [5]; alternatively, (E)^'^ = Dt, where D is the coefficient (20) in Ref. [5]. 



The fluctuating part (Ey t ' ' is given by Eq. (21) with [t/m] replaced by F t (m) = 
[t/m] — t/m. Clearly, (E)^ 1 '^ = identically only if F t (m) = 0, i.e., m — 1, for all relevant 
pairs (m, m'). This will occur in two cases: (a) The potential is mono-harmonic, i.e., V m = 
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for \m\ > 1. (b) The correlations Cp{m — m!) = for m ^ m', corresponding to a uniform 
(x-independent) distribution (13). In general, using |F t (m)| < 1 and (22), we easily see that 
< B = (k 2 C /2) J2mm>>o mm ' lKnKn'1; where C is an upper bound of \Cp{m)\ which 
decays with \m\. The quantity B is always finite for a typical, differentiable V(x) with \V m \ 
decaying faster than |m| -2 . Thus, (E)^' F ^ is bounded and is therefore much smaller than 
(E) t for large t. An example of significant fluctuations for small t is considered below. 



Exact closed results for a class of wave packets. We now focus on the class of wave packets 
for which ^(p) is piecewise constant on the unit intervals n < p < n + 1 for all integers 
n. This corresponds, by Eq. (4), to a /^-independent i/)p(x), i.e., the same i])p{x) = ip(x) is 
associated with all the (3-KRs. For /^-independent correlations in (13), Cp{m) = C(m), one 
has Cj(m) = C{m)5j$ in (14). Using the last result in the derivation above of Eq. (21), we 
see that Eq. (21) can be replaced by an exact formula for all times t: 



2 



(£>l 1) = y E mm'V m V^C(m-m') 

mm'>0 



t 



Also, in Eq. (18), Gp(m) = G (m) + (3C(m) for < (3 < 1. Then 



(23) 



Gj(m)= f d(3exp(-2irij(3)Gp(m) = 

J 



n , \ C(m) 
G (m) + 



*i,o + (1 - **o) • (24) 



2 

Inserting Cj(m) = C(m)Sj j0 and (24) in (19), we obtain 

h. t ( i \ lms 

(E)? ] = -£yE V m C(m) £ LJJ_. (25) 



2^ T^i s 

Let us briefly analyze the exact results (23) and (25). It is clear from the definition 

of m that the fluctuating part (E)^ 1 '^ of (23) is a periodic or quasiperiodic function of 
t depending on whether the number of harmonics of the potential is finite or infinite. 
Concerning (25), we first introduce the notations W\ = —k/(2irl)J2 m V2m+iG(2m + 1), 
W 2 = -k/(2irl) E m ^o^ m C(2m), and W = W l + W 2 . A good approximation of (25) for 
even I is then (E)f^ m W\h\{t) + 7], where 7 m 0.5772 is the Euler constant. For odd I, 
{E)f^ ks W^[ln(t) + 7] — Wi ln(2). Thus, (25) provides logarithmic corrections to the domi- 
nant growth of (E) t given by (23). These corrections vanish if, e.g., the symmetry conditions 



mentioned after Eq. (19) are satisfied. Then, (E) t = {E) Q + (E)\ so that the asymptotic 
behavior (23) of (E) t starts at t = 1. If, in addition, (E)f' F) = ( see the cases in which 
this occurs at the end of the previous section), (E) t is an exactly linear function of t. See 
also note [8]. 



In general, (Ey t ' may be significant relative to (E)\ ' for small t. As a simple example, 
we consider the case of the two-harmonic potential V(x) = cos(x) + r\ cos(2x) with 7pp(x) = 
tp(x) = A[l + Acos(x)], where rj and A are constants and A = [2n(l + A 2 /2)] -1 / 2 is the 
normalization factor, J 27r dx\ip{x)\ 2 = 1. The only nonzero values of V m and C(m) are those 
for \m\ < 2 and can be easily calculated from (6) and (13). Inserting them in (23) and in 
the results above for (25), we get 

(£>(« = e 

where F t (2) = or —1/2 for t even or odd, respectively, and (27) holds for even I (the 
case of odd I can be treated similarly and will not be discussed here). Choose A = V% 
so that |A|/(2 + A 2 ) assumes its maximal value a/2/4. Then, for i] 1 and sufficiently 
small k 1, the logarithmic correction (27) is significantly larger than (26) for all times t 
satisfying \n(t)/t 3> 2nlk. If k > 1, on the other hand, (E)f^ is negligible relative to (E)^. 
For T) ~ 1, the fluctuations in (26) have a magnitude of about 40% of that of the coefficient 
D ^ {E) ( p jt and may be therefore observed experimentally for not too large t. 



2 + A 2 



t + k' 



2r]X 
2 + A 2 



(26) 



Incoherent mixtures. We now assume that the initial KP state is an incoherent mixture 
of plane waves exp(ipx) with momentum distribution f{p) sufficiently localized in p. At 
time t, the expectation value of the kinetic energy in the state evolving from a plane wave 
with p = n + (3 is, from Ref. [5] (with a change in notation), 

E t (n,P) 



^ + ^E- 2 |^| 2 ^^. (28) 
2 r^o smW2) 
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The MKE of the mixture is then (E t ) = Jo^EJ( n + P) E t(n,(3). Defining the func- 
tion fo(/3) = J2 n f( n + 0)i which is periodic in j3 with Fourier expansion fo(/3) = 
J2j f(j) exp(27ri < 7/3), and using (15), we obtain after a straightforward calculation 

(Ek) = (E ) + k 2 Yl m 2 \V m \ 2 £ (-l) mla f(mla)(t - \a\). (29) 

m>0 a=l— t 

Clearly, in the asymptotic time regime t >> 1, (29) always exhibits a linear growth without 
fluctuations, unlike (21). For a nonuniform mixture [fo(/3) ^ 1, i.e., f(j) ^ Sj t0 ], this linear 
growth will be generally preceded by a transient behavior for sufficiently small t. 

In conclusion, we have shown that for KP wave packets the asymptotic quantum-resonant 
evolution of the MKE is typically given by (E) t ~ Dt + (E)f ,F \ where D is the coefficient 
of linear growth [5] and (E)f ,F ^ are bounded fluctuations of a simple number-theoretical 
nature. For a large class of initial wave packets, (E) t is exactly the superposition of its 
asymptotic behavior and transient logarithmic corrections which may be negligible for large 
k or can vanish under some conditions. The asymptotic behavior then starts already at small 
t, where the fluctuations are most prominent, especially when the magnitude of (E)f ,F ^ is 
comparable to that of D as in the example above. We remark here that the sensitive depen- 
dence of D on the harmonics of V(x) and on the initial wave packet, pointed out in Ref. [5], 
is featured also by (E)\ 1,F \ as one can clearly see from (21) or (26). The fluctuations in (E) t 
for small t may be noticeable even when the transient corrections are significant. All the 
quantum-resonance behaviors for t not too large are expected to be robust under sufficiently 
small deviations of r/(2n) from integers and may be therefore observed experimentally. 

Since integer values of t/(2tt) correspond to strong quantum regimes, the fluctuations 
and transient phenomena studied here are basically different in nature from known ones 
occurring in semiclassical regimes [9-11], e.g., the asymptotic quasiperiodic fluctuations 
associated with dynamical localization and transient behaviors featuring an approximate 
coincidence of the classical and quantum evolutions on some initial time interval. It would 
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be interesting, however, if one could establish some correspondence between the general 
quantum-resonance phenomena in this paper and classical behaviors of related systems, as 
it was done in a particular case in Ref. [4]. 
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